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While electrowetting has many applications, it is limited at large voltages by contact angle sat- 
uration — a phenomenon that is still not well understood. We propose a generalized approach for 
electrowetting that, among other results, can shed new light on contact angle saturation. The model 
assumes the existence of a minimum (with respect to the contact angle) in the electric energy and 
accounts for a quadratic voltage dependence ~ in the low-voltage limit, compatible with the 
Young-Lippmann formula, and a ~ saturation at the high- voltage limit. Another prediction is 
the surprising possibility of a reversed electrowetting regime, in which the contact angle increases 
with applied voltage. By explicitly taking into account the effect of the counter-electrode, our model 
is shown to be applicable to several AC and DC experimental electrowetting-on-dielectric (EWOD) 
setups. Several features seen in experiments compare favorably with our results. Furthermore, the 
AC frequency dependence of EWOD agrees quantitatively with our predictions. Our numerical re- 
sults are complemented with simple analytical expressions for the saturation angle in two practical 
limits. 



I. INTRODUCTION 

The term Electrowetting, in its broadest sense, refers to 
techniques by which one can control the apparent wet- 
tability (characterized by the contact angle) of liquids, 
by applying an external electricpotential While 
it has numerous applications |12| - [20| . electrowetting is 
known to be limited by the so-called contact angle satu- 
ration (CAS) [i,0,|2lHl3 as depicted in figure [H As the 
term indicates, an electric potential can incur a change in 
the contact angle, but only to a certain extent. Further 
voltage increase has no additional effect on the contact 
angle. This behavior is not accounted for by the stan- 
dard model of electrowetting, and its origin still remains 
a point of controversy 0-0) fisl . [2l| - [33j . 

When a small drop of liquid is placed on top of a solid 
surface it assumes the shape of a spherical cap [343 • The 
contact angle 6*0 between the drop and the surface, given 
by the Young formula [33 - [36| 
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depends on the three interfacial tensions: solid/air 7sa, 
solid/liquid 7si and liquid/air 71a, where the air phase 
can be replaced by another immiscible fluid [tI. IT7l.1l9j. 

The Young formula, eq [U can be obtained by mini- 
mizing the capillary free energy with a fixed volume con- 
straint 

i^cap (f ) = ^a7sa + ^l7sl + ^la7la - V AP (2) 
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FIG. 1. A schematic plot of electrowetting contact angle, 0(U). 
For zero voltage, the contact angle is the same as the Young 
angle, ^o- At low applied voltages U , the contact angle follows 
the Young-Lippmann formula, eq [6] (dashed line) cos6{U) — 
cos6o — cos9yl{U) — cos^o ~ , but for higher voltages 
the contact angle gradually deviates from the Young-Lippmann 
behavior and saturates towards some finite value, Ssat- 



where Aij are the interface areas between the i and j 
phases, z, j = a (air), 1 (liquid), and s (solid); the drop vol- 
ume is V and the pressure difference across the liquid/air 
interface is AP. For partial wetting, 7sa < 7si + 7iaj the 
capillary free energy Fcap has a minimum at the Young 
angle, 6*0 [figure [Ifa)] . 

The contact angle 9 can be varied from its initial 
value 9o by applying an external voltage of several volts 
to several hundreds of volts across the liquid drop. 
A commonly used electrowetting setup developed by 
Rinkel et al [l[ and later perfected by Vallet et al Q 
is called electrowetting-on-dielectric (EWOD). The ap- 
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FIG. 2. Schematic plots of (a) the capillary free energy _Fcap for a 
spherical drop of fixed volume and (b) a hypothetical (negative) 
total capacitance, —Ctot, as function of contact angle 6 in elec- 
trowetting setup. The minimum of the capillary term Fcap occurs 
at the Young angle 6o, while the electric term F^i — —^CtotW^ 
(or equivalently of —Ctot) is assumed to have a minimum at a 
finite saturation angle, Ssat. 



paratus, roughly sketched in figure [3{a), includes a flat 
electrode as a substrate, which is coated with a thin di- 
electric layer (tens of nanometers to several micrometers 
thick), whose purpose is to prevent Faradic charge ex- 
change {i.e., electrochemical reactions) at the electrode. 
It is common that this dielectric layer is then topped with 
an even thinner hydrophobic (i.e., Teflon) layer in order 
to control its surface tension. A drop of ionic solution 
is placed atop the coated electrode and a thin counter- 
electrode (usually a bare platinum fiber) is inserted into 
the drop from above. The drop is surrounded by air or by 
another immiscible dielectric liquid. Applying a voltage 
across the drop can cause a large change, of several tens 
of degrees, in the contact angle. 

As reviewed in Refill a simple relation between the 
contact angle and the applied voltage can be derived. 
When an external voltage U is applied an electric double- 
layer is formed at the liquid/substrate interface. The 
total free energy Ftot has two contributions: a capillary 
term i^cap, defined [131 in eq[21 and an electric term, F^i 
that depends on 9, U and other system parameters. 

Ftot(0,f/) -Fcap(e)+Fel(0,C/) (3) 

Within the standard model of electrowetting (under ex- 
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FIG. 3. (a) A sketch of an EWOD setup with AC voltage U of 
frequency /, and (b) its equivalent AC circuit. The two parallel- 
plate capacitors of capacitances Ci and C2 and areas Ai and A2 
represent the substrate/liquid and the counter-electrode/liquid 
interfaces, respectively. The bulk liquid drop is represented by 
a resistor Ri through which the two capacitors are charged and 
discharged. The charging time through the resistor is equal to 
the build-up time rt of the double-layers at the two interfaces. 
The discharge resistor R2 represents the Faradic charge transfer 
processes that relax the electric double-layer near the counter- 
electrode with relaxation time Tr. Such mechanism is prevented 
at the substrate electrode because of its dielectric coating. 



ternal voltage control) the electric term is evaluated as 

1 



(4) 



where Cm is the capacitance of the liquid/substrate in- 
terface, which is modeled as a parallel-plate capacitor: 



Cid — £o^u 



£1 



£o£d 



(5) 



where is the substrate area that is covered by the 
liquid drop, d is the width of the dielectric layer and Sd 
and El are the dielectric constants of the dielectric coat- 
ing and liquid drop, respectively. The Debye screening 
length, Ad, is the width of the electric double-layer. In 
most cases, d/ed ^ Xd/ei, and the second approxima- 
tion in eq \E\ can be justified. The parallel-plate capac- 
itor model has been shown to be valid as long as 
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fringe fields are negligible {d ^ A^^ ) and the drop of 
volume V is not too small for double-layers to be cre- 
ated (Aj) <^ V^^'^). An implicit assumption made in eq|3] 
is that the counter-electrode does not contribute to the 
total capacitance. 

In an electrowetting setup, the contact angle 0{U) de- 
pends on the external voltage U. By substituting eq [5] 
into eq H] and minimizing i^tot of eq [3] with a fixed vol- 
ume constraint, the Young-Lippmann formula for the 
contact angle Oyl{U) is obtained 

cosOYhiU) = cos9o+g-^U^ (6) 

where = eoed/i^-find). It is a common practice to 
extend this DC voltage model to AC setups by using the 
rms voltage in eqlHl U^ms- Note that similar results 

can be obtained using force balance at the three-phase 
contact line [38| . 

Experiments [1-0, [H, HH-lll have shown that 

the ~ behavior predicted by the Young-Lippmann for- 
mula is indeed found for a range of low applied voltages, 
but the pre-factor of the U'^ term, eq[Bl does not usually 
match the experimental data. For larger values of U, a 
deviation from the behavior is observed and a satura- 
tion in the contact angle, 9{U) -> ^sat, is reached gradu- 
ally, as is schematically sketched in figured] In addition, 
it is convenient to define a characteristic value of the 
cross-over voltage U* by requiring that 9yl{U*) = 9sat 
in eq[6] 

([/*)' =g(cos0,at- cos 0o) (7) 

Over the last decades, several models have been pre- 
sented in an attem pt t o explain CAS jl^] . Most of these 
models 0, [H, [H, Isil - lssj are based on specific leakage 
mechanisms. Others, as in Ref [13, proposed heuristic 
arguments in order to predict CAS in electrowetting sys- 
tems without relying on a specific mechanism. 

Considering that the origin of CAS is not well under- 
stood from general principles, the objective of the present 
work is to offer a different approach to CAS, and to elec- 
trowetting in general. In the following section we consider 
the general circumstances in which CAS can occur intrin- 
sically (without leakage). We present a low- voltage limit 
compatible with the Young-Lippmann quadratic voltage 
dependence and a high- voltage limit in which CAS is ob- 
tained. Furthermore, we identify a possibility for a novel 
electrowetting regime we call reversed electrowetting. In 
section III we present an application of this approach to 
EWOD experimental setups using a geometry-dependent 
model, and use AC circuit analysis to calculate the free 
energy. Section IV is dedicated to showing several numer- 
ical and analytical results and their compatibility with 
experiments. We conclude in section V with a summary 
along with some further discussion and an outlook on 
future research. 



II. A GENERALIZED MODEL OF 
ELECTROWETTING 

A. Generalized free energy and contact angle 
saturation 

Our starting point is eq [3] above. Assuming that 
all the electric energy is stored via charge separation, 
it can be written in terms of the total capacitance, 
Foi — —\Cto'JJ^- The total electrocapillary free energy 
Ftot is now written as 



Ftot{9, U) = Feap(0) - ^Ctot{9)U^ (8) 

where all capillary contributions are included in i^cap- 
As Fcap is independent of U , the relative magnitude of 
the two terms in eq |5] is controlled by the negative 
dependence of the F^y term. 

Physical insight can be gained from eq [8] by making 
different assumptions regarding the behavior of Ctot(^)- 
Particularly, interesting results follow from the assump- 
tion that Ctot(^) has a maximum at some finite angle 
(see figure [Sfb)), as will be shown below to be the case 
for EWOD setups. Our model therefore differs from pre- 
vious models that took Ctot to be equal to Cm, (eqU]), 
which yields a monotonically decreasing Fc\. For reasons 
to be immediately apparent, we denote the angle where 
C'tot(^) has a maximum (or, equivalently, F^x has a min- 
imum) as Osat- This angle, in general, is different from 
the Young angle, 9q, which minimizes the capillary term, 

^ cap • 

We now show how the existence of a global electric free- 
energy minimum at a finite contact angle yields CAS. 
With no applied voltage {U ~ 0), i^tot — ^cap and 
the system adheres to the Young angle, 9 = 9o, (fig- 
ure [DJa)). Similarly, when the applied voltage U is very 
large, the free energy is dominated by the electric term, 
I Foil K Fcap, and the system tends towards ^sat, 

which minimizes F^i (or, equivalently, maximizes Ctot, 
figure IDJb)). Now, if the two contributions are concave 
for the accessible range of 9, then the minimum of Ftot 
shifts smoothly from 6*0 towards 6'sat as U is increased 
from zero to an arbitrary large value, as is schematically 
illustrated in figure 21 This description is consistent with 
CAS and implies that the saturation angle found in ex- 
periments can be identified with our definition of Sgat • 

Below we analyze more quantitatively the conse- 
quences of such a global electric minimum at the low- 
and high-voltage limits. In the former, we show a ^ C/^ 
variation of the contact angle with a pre-factor that can 
match the Young-Lippmann formula or be different from 
it. In the latter, an asymptotic ~ C/~^ approach to ^sat 
is found. 
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FIG. 4. A schematic plot of the total free energy Ftot{0,U) = 
Fca,p{9) + Fi,i{9, U), for a sequence of five applied voltages: U4 > 
Us > U2 > Ui > U = 0. Since the two energy terms are 
concave in the accessible range of 6, their sum is also concave 
and has a minimum. Increasing U from zero to a large f/4 causes 
a gradual shift of the minimum from 9o towards Ssat. Note that 
it cannot cross beyond 6sat regardless of how high the voltage U 
is because 6sat is the minimum of the dominating Fd term. 



change its sign (see section II. D). It is shown in section 
IV. D under which conditions the pre-factor converges to 
that of the Young-Lippmann formula for low voltages in 
typical EWOD experimental setups. 



C. The High Voltage Limit 

For [/ ^ 00, the electric energy becomes large relative 
to the capillary energy and so the minimum of Ftot occurs 
at e{U) = 9sat + S9. Expanding the condition Fl^^{0) = 
around flgat, one has 



F' I 

cap 
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e{u) 



F' ( 

cap \ 



(13) 
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Ctot(^sat) 



(14) 



Hence, saturation in 6 is approached asymptoticall y, a s 
in qualitative agreement with experiments Il4|. 



D. Reversed Electrowetting 



B. The Low Voltage Limit 

For J7 0, the minimum of i^tot occurs close to 6q. 
Expanding this minimum condition Fl^^^{9) = to first 
order in 59 = 6{U)~9o, while recalling that i^cap(^o) = 0, 
we obtain 

^^c'ap (^0) 69 - hj^ \C[,, (0o) + C'U i9o) 69] ^ (9) 
yielding, 

''-2F^'^^i9o)-^,U^Cl'^A0o)'' ^''^ 
and to leading order in one has 

9iU)^9, + l^if^u^ (11) 
and equivalently 

cos^(L/).cos^o-i%|^C/^ (12) 

^ ^cap 1^0) 

We see that at low voltages, the deviation from the 
Young angle is proportional to C/^, just as in the Young- 
Lippmann formula. However, the pre-factor is a function 
of ^0 and can take different values than in eqlHl and even 



An interesting conclusion can be drawn from the dis- 
cussion in section II. A. Recalling that in our model elec- 
trowetting results from an interplay between capillary 
and electric energies (each with its own minimum at ^0 
and 0sat, respectively), as voltage is increased the electric 
energy gradually becomes dominant and the contact an- 
gle is driven away from flp towards flgat (figure |4]) . Since 
^0 is determined only by the capillary parameters (as in 
the Young formula, eq[T|) and 6'sat is determined solely by 
the electric parameters, it is possible to envisage a sys- 
tem in which the saturation angle flgat is actually larger 
than the Young angle, 6'sat > ^'o rather than 6'sat < ^'o as 
in the usual case. In such a setup, applying a voltage will 
cause an increase of the contact angle, in total contra- 
diction with the Young-Lippmann formula, eq[6l Hence, 
the model proposed here allows for the possible existence 
of a new regime of electrowetting, which we refer to as 
reversed electrowetting. 

By examining the slopes of each energy term near the 
minimum of the other (see figure [5]), it is possible to show 
that in the low- and high-voltage limits, eqs [TT] and [TH 
the pre-factors of both and terms can take ei- 
ther positive or negative values depending on whether 
9q > 6'sat or 9o < 6'sat; for the low- voltage limit F^ap (6'o) 
is positive by definition, but C'^^^ {9o) is positive only if 
9o > 9 sat and negative for < ^sat- Likewise, for the 
high- voltage limit (7"^^ {9 sat) is negative by definition but 
Fp'ap {9sat) is negative only if 6'o > 0sat and positive for 
^0 < ^sat- Thus, we have shown how reversed electrowet- 
ting manifests itself in those limits. 
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III. A TWO ELECTRODE MODEL OF EWOD 

Our goal in the remainder of this work is to elaborate 
on the physical conditions that are involved in determin- 
ing a finite ^sat angle in specific EWOD experimental 
setups. However, we would like to stress that the pro- 
posed mechanism is general and may be applied to other 
realizations and experimental setups manifesting CAS. 



A. System Setup &z Geometry 

A setup of an EWOD setup is presented in figure [5] 
The drop (dielectric constant £;) is assumed to retain its 
spherical-cap shape, with height h from the surface, to- 
tal volume V and contact angle 9. The metal electrode 
is coated with a dielectric layer of thickness d and dielec- 
tric constant Ed- The top counter-electrode is modeled 
as a thin cylinder of radius b and the gap between the 
two electrodes is hg. The two electrode areas covered by 
the liquid are Ai and A2, respectively. For spherical-cap 
shaped drops Ai and A2 are related to the contact angle 
6 through the fixed volume constraint: 



Ai 






A2 








tan — 
2 


h 

a 




V 


= ^(3a 


+ h^) 



(15) 



where a is the radius of the covered portion of the sub- 
strate electrode. 

It should be noted that can only take values in the 
range 6*111111 < 9 < it. The lower limit, 6*111111, occurs when 
the drop height matches the gap between the lower tip 
of the counter-electrode and the substrate, h — hg (see 
figure [5|), resulting in 



cot ■ 



/ 2V 



(16) 



We will show that 9s, 



> 



hence 6'„iin is an inacces- 



sible lower bound of the contact angle. The upper limit 
9 = TT is the de- wetting limit. 



B. The AC Free-Energy 

The free energy, eq[51 depends on the total capacitance 
Ctot, which includes all relevant contributions. Unlike 
the traditional Young-Lippmann treatment in which only 
the capacitance of the liquid/substrate interface is taken 
into account, we consider explicitly the existence of an 
additional double-layer, residing at the interface between 
the liquid drop and the counter-electrode. 

Experimental setups and applications usually employ 
AC circuits to produce an electrowetting effect. Under 




FIG. 5. A schematic EWOD setup as used in our analysis. A 
liquid drop shaped as a spherical cap of volume V, height h and 
dielectric constant ei is placed atop a flat metal electrode. The 
metal electrode is covered with a dielectric coating of thickness 
d and dielectric constant Sd- A metal wire (used as a counter- 
electrode), modeled as a thin cylinder of radius b, is inserted into 
the drop from above. The gap between the two electrodes is hg. 
The area of the substrate electrode covered by the drop is Ai 
and that of the counter-electrode is A2. The applied voltage is 
U and the contact angle with the substrate is 6. 



those circumstances, double-layers are transient: with 
each AC half-cycle a double-layer of opposite polarity 
is formed at each electrode/liquid interface and subse- 
quently dissolved away. In addition to relaxation by re- 
versal of polarity, other mechanisms of relaxation can 
act at the counter-electrode/liquid interface such as elec- 
trochemical Faradic processes. The dynamical processes 
are, therefore, governed by two intrinsic time scales (be- 
side the AC frequency): 

1. The double-layer build-up time, Tb, which can be 
estimated to be 



n 



D 



(17) 



where A_d = y/ eQeiksT /2e^Csa,it is the Debye 
length, ksT the thermal energy, Cgait the salt con- 
centration, D is the diffusion constant and _L is a 
typical system size 39]. 

2. The double-layer relaxation time, r^, which can 
similarly be expressed in terms of system param- 
eters through the RC circuit relaxation formula 



Tr = R2C2 = 

A-2 



(18) 



where p — R2A2 is defined as the zero-current 
(Faradic) resistivity to charge transfer by electro- 
chemical processes and A is the contact area. 
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In order to discuss the period-averaged properties of 
the system, we employ a standard AC circuit analysis. As 
shown in figure ISjb), we model the two liquid/electrode 
interfaces as two capacitors with capacitances Ci, C2, de- 
fined in a similar fashion as in eq[5] 



C2 



A 



D 



(19) 



Note that the main contribution to Ci comes from the 
coated dielectric layer of thickness d {d/si :s> )^d/£i), 
while for C2, the only contribution comes from the double 
layer of thickness Ad (because the counter-electrode is 
not coated). The cylindrical geometry of the counter- 
electrode is not considered because & 3> A^. 

The two capacitors are charged and discharged through 
a resistor Ri that represents the bulk of the liquid 
drop. The relaxation of the double-layer at the counter- 
electrode is modeled by an extra discharge circuit with 
a resistor R2, while the capacitor Ci does not have a 
discharge circuit since charge transfer at the substrate 
electrode is prevented by its dielectric coating. The ap- 
propriate resistance values can be inferred from the build- 
up (rf,) and relaxation (r^) times, egs flTl and [TSl again 
through the RC circuit relaxation formula: 



Ri = rb(Cr 
R - 



(20) 



Drawing on the AC circuit analogy, the period- 
averaged free energy is: 



Ftot(^, U, U) = i^cap + Fel = i^cap(0) 



1 



1 



2uj Ztot(0,uj) 



f 

(21) 

where U is understood to be the rms value and Zt^t is 
the total impedance of the circuit, [figure Eljb)], which 
can be represented schematically as 



Ztnt — Zr 



Zr,(S{Zc,\\Zr,) 



(22) 



It is straight forward to show that the squared magnitude 
of the total impedance is 



1 
1 



CIC2 

-2ny\ 



TrjTr + Tfc) 
1 + w2r2 



(23) 



Since Ci is proportional to Ai , which vanishes at — > tt 
and C2 is proportional to A2, which vanishes at 6 ^ ^min, 
\Ztot\ diverges both at 6* — 6'min and 6 tt. Therefore, 
it must have a minimum at some intermediate value: 
< ^sat < TT. Heucc, our model of electrowetting 
presented in section II. A is indeed applicable to typical 
EWOD setups. 

Substituting eg s [T5] and [T^ into egs and vields an 
expression for Ci (9) and C2 (0) and, consequently, for Fc\ 



as a function of 9. Its minimization can be done numer- 
ically (section IV. A) and yields the equilibrium contact 
angle (for given applied voltage and frequency) . In some 
limits (sections IV. C and IV. E) analytical approxima- 
tions can be derived as well. 



IV. RESULTS AND DISCUSSION 

A. The Electrowetting Curve, e{U) 

In order to demonstrate quantitatively the model va- 
lidity, we performed numerical calculations for parameter 
values that are in accord with some typical experimental 
setups. In figure [6] we present the reactance 1/ujZ (ap- 
pearing in eq I21[) computed for parameter values as in 
table H] and with an AC frequency / = u!/2tt = IkHz. 
The build-up time was calculated using eq [T7] to be 
Tf, — 1.34 ms, while the relaxation time was calculated us- 
ing eq[TH]with p = 1 17 • m^, yielding [H] = 0.53s. 
For the chosen values of parameters, the ratio of ca- 
pacitances for zero voltages {6{U = 0) = ^o) is about, 
C2/C1 ~ 25. In the figure a maximum at a finite angle 
Osat = 53.3° is clearly seen. Notably, this saturation an- 
gle is much larger than the minimal possible angle in this 
setup, ^min = 37° (scc cq [T6|) . As a consequence CAS is 
obtained for finite values of A2, much before the limit 
A2 ^ characteristic to ^min- 
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FIG. 6. Total reactance l/(aj|^tot |) (in nanoFarad) as a function 
of the contact angle. Parameters for the EWOD system are 
Tt — 1.34 ms, Tr — 0.53 s and / = ijj/2it = 1 kHz, and all 
other parameters are chosen as in table [T] The reactance has a 
maximum for a finite value of 6sat ~ 53.3°, at which the system 
exhibits CAS. The minimal contact angle is Omin = 37°, eg 1161 



Figure [7] presents the calculated electrowetting curve 
0{U) for the same system, where 0{U) is calculated by 
minimizing i^tot from eq 1211 together with a plot of 
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Parameter 


symbol 


value 


dielectric constant of liquid 


£l 


80 


Debye length in liquid 


Xd 


1.34nm 


volume 


V 


5/iL 


width of dielectric layer 


d 


0.1/im 


dielectric constant of dielectric layer 


£d 


2.67 


liquid/air surface tension 


7la 


72.8mN/m 


dielectric/air surface tension 


7sa 


12.7mN/m 


liquid/dielectric interfacial tension 


7sl 


47mN/m 


gap between counter-electrode and substrate 




0.7mm 


radius of counter-electrode 


b 


12.5/^m 



TABLE I. Parameter values of a typical elect rowetting setup. 
The liquid drop contains an aqueous ionic solution and is placed 
on top of a Miyaline-C/Teflon substrate. 
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the Young-Lippmann formula (dashes), where an effec- 
tive g'^^ pre-factor is used to fit the full calculation. 
This is similar to what is done in many experimental 
works where the g value is fitted from the low U de- 
pendence, and not by using explicitly eq [6l We use a 
specific = (1 + uj^T^)~^^^ as derived in section IV. D. 
The figure shows that several common experimental fea- 
tures are reproduced (as compared with the schematic 
figure [1]); an initial compliance with the scaled Young- 
Lippmann formula at low voltages is followed by a cross- 
over at intermediate voltages to a different regime. Using 
eq □ (with /ff), the cross-over voltage is evaluated to be 
U* ~ 76.9 V. At U > U* an asymptotic convergence 
of the contact angle towards a saturation value is seen, 
0{U) 9sat = 53.3°. This is further demonstrated in fig- 
ure [71 where the asymptotic 9{U) —9sat ~ U~'^ is plotted 
(dotted line) following eq [TH The asymptotic behavior 
approximates rather well 9{U) for voltages larger than 
120V. 

It is appropriate to define another voltage, [/gat, char- 
acterizing the saturation range of the potential. An op- 
erational definition that we employ is that at L'sat, the 
calculated 9{Usa.t) deviates from 6'sat by 2%. With this 
definition, we obtain Usat — 252.1V. The electrowetting 



curve presented in figurej 
perimental observations 



_agrees qualitatively with ex- 
13 . [28j , which show the effect 
of contact angle saturation. Unfortunately, because the 
parameter values needed for quantitatively comparison 
with experiments lack at present, we used instead rea- 
sonable estimations. 



B. The Reversed Electrowetting Curve, 6{U) 

We now illustrate how reversed electrowetting ^sat < 
^0, which is a natural outcome of our model, can be seen 
in the laboratory. Let us consider a system similar to 
the one presented in the previous section with the two 
following changes (see table |lTl : the interfacial tensions 



FIG. 7. Calculated contact angle as a function of applied volt- 
age, U (full line), as well as a plot of the Young-Lippmann for- 
mula using g'^^ (dashed line, see text). Parameter values of 
the EWOD system are taken from table [T] and rt — 1.34 ms, 
Tr = 0.53 s and / = u)/2-k = 1 kHz. Many of the features 
of CAS, shown schematically in figure [l] are reproduced. At 
U = Q, the contact angle is 6q = 118°. For small U , there is an 
initial compliance with the rescaled Young-Lippmann formula, 
eq [6] 0{U) ~ OYh{U), followed by a cross-over occurring at 
U* — 76.9 V, calculated from eq[7] At larger U, the contact an- 
gle tends asymptotically to a saturation angle ^sat = 53.3°. The 
asymptotic 6{U) — Ssat ~ is plotted (dotted line) following 
eq[l4]and approximates rather well 9{U) for voltages larger than 
120 V. An operational definition of the saturation voltage (/sat 
(see text) yields Usat — 252.1V. Note that Osmt is conceivably 
larger than the minimal possible angle Omin = 37°, eg 1161 



are chosen such that — 60° , and we now model a non- 
polar liquid with a dielectric constant e; = 2, which yields 
a build-up time of Tf, = 1.34 ms and a relaxation time of 
Tr — 13.2 ms. The AC frequency is / = a;/27r = IkHz as 
before. 

Figure [5] presents the calculated electrowetting curve 
9{U). The plot features an initial compliance with the 
negatively rescaled Young-Lippmann formula (at low 
voltages such that the contact angle increases with the 
applied voltage. This is followed by a cross-over at 
intermediate voltages towards saturation. Using eq [7| 
(with g"^), the cross-over voltage is evaluated to be 
U* ~ 65.5 V. For U > U* an asymptotic convergence 
of the contact angle towards a saturation value is seen, 
9{U) 6'sat = 101.9°. Using the same definition as 
in section IV. A, the saturation voltage is found to be 
f/sat ^ 208.5 V. 

Since our reversed electrowetting predictions (figure [7]) 
are rather for specific parameter values, it will be of ben- 
efit to check their validity with experiments conducted 
on similar electrowetting setups. 
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Parameter 


Value 


ei 


2 


\d 


1.34 nm 


V 


5 i-iL 


d 


0.1 /im 


Ed 


2.67 


7la 


20 mN/m 


7sa 


15 mN/m 


7sl 


5 mN/m 


hg 


0.7 mm 


b 


12.5 fim 



TABLE 11. Parameter values of a hypothetical system that ex- 
hibits reversed elect rowetting. 




300 



ration angle reaches a constant value for the entire high 
frequency range down to / ~ 1 kHz, even for the smallest 
of the chosen relaxation times (r^ = 5.3 ms). Moreover, 
the larger is, the wider is the range for which the 
saturation angle is constant. This can be explained by 
taking into account that whenever ^ 1//, the counter- 
electrode double-layer can hardly relax. 




10'^ 10'^ 10° 10^ 10^ 10^ 10'* 



U[V] 



f[Hz] 



FIG. 9. Saturation angle, Osat, as function of the AC frequency 
/, for the EWOD system with parameter values as in table [T] 
Tt = 1.34ms and for several values: Tr = 0.53s (solid 
line), 0.053s (dashed line), 5.3 ms (dashed-dotted line). In the 
range / > 1 kHz, the AC saturation angle, which takes into ac- 
count relaxation of the double-layer by electrochemical processes, 
matches the r,. — > oo (dotted line) of — 53.3° even for the 
smallest r,. = 5.3 ms (fastest Faradic relaxation). 



FIG. 8. Calculated contact angle as a function of applied volt- 
age, U, (solid line) for reversed electrowetting together with a 
manually scaled {g°^ < 0) Young-Lippmann formula (dashed 
line). Parameter values of a EWOD system are taken from ta- 
ble [TO and Tb — 1.34 ms, Tr — 13.2 ms and / — uj/2-k = 1 kHz. 
At U — 0, the contact angle is 9o = 60°. Even for small U, 
there is a deviation from the naive Young-Lippmann formula, 
eq[Sl because 6{U) — 6o ~ U^, with a positive pre-factor. The 
rise is followed by a cross-over occurring at U* — 65.5 V. At 
larger U, the contact angle tends asymptotically to a saturation 
angle ^sat = 101.9°. Using the same definition of the saturation 
voltage (7sat as in figure|7]yields t/sat — 208.5V. 



C. The Frequency Dependence of Electrowetting 

In order to explore the effect of the frequency of the AC 
voltage within our model, the dependence of the satura- 
tion angle ^sat on frequency was calculated numerically 
by minimizing eq [5T] and plotted in figure IH] for several 
values of = 0.53 s, 0.053 s, and 5.3 ms. It can be seen 
that for this specific choice of parameters, the AC satu- 



The frequency dependence of the; contact angle has 
been experimentally studied in Ref 8 for several applied 
voltages. Figure [10] shows a comparison of a minimiza- 
tion of eq[23]for a range of frequencies, with experimental 
data. The calculations have been performed for system 
parameters as in table IIIIl which have been inferred [45| 
from Ref 1^, except for A/^ = 300 nm, which was used 
as a fitting parameter to the experimental results. This 
value corresponds to ionic strength of less than 10~^ M 
and is compatible with de- ionized water used in Ref |^. 
The build-up time was deduced from the featured experi- 
mental results to be ti, = 0.1 ms, and the relaxation time, 
Tr = 2.4ms, was calculated using eq[T5]with p = 1 17 -m^. 

Comparison of the two plots shows that our model re- 
produces rather well the frequency dependence found in 
experiments. It can be seen that the electrowetting effect 
diminishes with rising frequency, and seems to vanish at 
/ > 100 kHz where there is hardly any deviation from the 
Young angle. This can be easily understood taking into 
account that for / > 10 kHz, 1/f becomes small as com- 
pared to the double-layer build-up time Tb — 0.1 ms. Un- 
der those circumstances, ions move too slowly and cannot 
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Parameter 


Value 


ei 


80 


\d 


300 nm (fitted) 


V 


5 /iL 


d 


5 /im 


Sd 


2.67 


7la 


72.8 mN/m 


7sa 


12.7 mN/m 


7sl 


47 mN/m 


hg 


0.7 mm 


b 


40 fim 



TABLE III. Parameter values of an electrowetting system of 
de-ionized water solution on Miyaline-C/Teflon substrate. The 
parameter values are inferred from Ref 8, except for Ad, which 
was fitted to obtain a quantitative agreement with experiments. 
The resulting value is indeed compatible with that of de-ionized 
water (ionic strength less than 10"" M) |'45^] . 



build considerable over-concentrations at the electrodes. 
Thus, at those high frequencies the period-averaged effect 
of the double-layers decreases considerably. 

The above results imply that it is of value to further 
explore the limit of slow electrochemical processes, ^ 
Tf,, uJTr 3> 1, which leads to a simplified expression for 
the free energy. Substituting eq [22] into eqUHwe obtain 
in this limit 



CtotiO) 
1 



2uj Zu 



2^1 



|:C/2(24) 



where = Cf ^ + is the total capacitance of an 
equivalent system without any electrochemical processes. 
The (1 -I- w^T^)"^/^ pre-factor in eq [24] depends on the 
AC frequency and reflects a diminishing electrowetting 
effect for rising frequencies. 

We note that the DC limit can be obtained by first 
assuming no electrochemical processes, — )• oo in eg [23l 
(leading to eq lM]) . and only then taking the DC limit of 
w to get Fei = —\CtotU^ of eq[5] This limit can 
be useful in applications where both the substrate and 
counter-electrode are dielectrically coated. 

If electrochemical processes are not totally excluded, 
but are just slow (order of seconds, in accordance with the 
value calculated in section IV. A), it is expected that the 
results obtained in this work will be applicable for that 
time scale, above which other mechanisms might take 
over. Such time-dependent behavior has been observed 
in Ref |2J. 
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FIG. 10. Contact angle as function of AC frequency / for several 
applied voltages U . (a) Calculated values corresponding to 57 V 
(solid line), 93V (short dashed line), 113 V (long dashed line) and 
143 V (dashed-dotted line) with parameters as in table fTTTI and 
Tb = O.lms, Tr — 2.4 ms. The value of 6o — 118° is indicated by 
a dotted line, (b) Experimental results adapted from Ref 8 for the 
same voltages as in (a): squares (57V), circles (93V), triangles 
(113 V) and inverted triangles (143 V). A quantitative agreement 
between the calculation and experiment can be seen. Note also 
that the electrowetting effect diminishes at high frequencies as 
predicted analytically in eq 124! 



D. Convergence to the Young-Lippmann Formula 



In its DC limit (w — >■ 0), equation [24] provides a path- 
way to estabhsh a relationship between our model and 
the Young-Lippmann formula and to show the conditions 
under which the two converge. Using eqs 1151 and 1191 (see 
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Appendix for more details), we have 



Acute saturation angles (large /3 ) 



Ctot — Cl 







hi 









-1 2-e 



where I = (SV^/tt)^/"^ is a typical drop length, 
is a dimensionless parameter, and 



1 — cos ( 



2 + COS 61 



1/3 



(25) 
(26) 

(27) 



is a monotonically increasing function of < < tt. 

As long as the second term in the brackets of eq[25]is 
small, our model agrees with the standard model, -Fci — 
— ^CiC/^, eqSJ with Ci = Cm- For a typical system, as 
the one presented in table I, the value of the constant pre- 
factor is rather small /3~^ ~ 0.01. Since hg/l < ^ < 23 , it 
is clear that the quotient can only be large when ^ hg/l 
or, equivalently, when the contact angle becomes small 
enough. Otherwise, the second term is negligible and the 
two models converge. 

Note that this view of the validity of the Young- 
Lippmann formula as being related to a certain range 
of the contact angles is a departure from the common 
approach which regards its validity being related to a 
certain range of applied voltages. 

By creating this link between the Young-Lippmann for- 
mula and our model it can be deduced that the proper 
way of extending the Young-Lippmann formula (within 
its validity range) from DC to AC is to replace — J> 
^mis/ + ^^'Tb ■ This is exactly the how the Young- 
Lippmann formula was scaled (by g'^^) in section IV. A. 



E. The Saturation Angle for Slow Relaxation 

{ujTr » 1, Tr » Tb) 

Within the slow relaxation framework, eg 1241 the mini- 
mization of -Fci (yielding 9sat) is equivalent to minimizing 
the total inverse capacitance a = 1/ Ctot • Using eqs 1151 
and fT9l we obtain 



a(0 



'Sd 



2-KeQedP 



1 



2/3- 



lf2 
9S 



i-i-^K 



(28) 



Minimizing a (^) yields a 6* order polynomial in ^ 

- 2 (2 + pr'^hl) f ~ 2/3^^ 



(29) 



It is possible to examine two separate limits for minimiz- 
ing a (^) , leading to two simple analytical expressions for 



If {Osat) ^ 2 then, near its minimum, eq [28] reduces 



to 



2/3- 



with a minimum at ^ = ^sat that satisfies 

-2 



I 



i-l-^hg 



= 



The solution yields 



and 0sat can now be obtained 



cos 6'sat = 



1 - 2dt 

1+4'at 



(30) 



(31) 



(32) 



(33) 



Inserting typical values from table U to check for self- 
consistency, we get ^sat ~ 0.173 <?C 2 as required. Using 
eq [221 the saturation angle in this case is calculated to be 
6'sat — 56.1°, which is not far from the value obtained by 
a full numerical calculation, 53.3°. As a rule of thumb 
we remark that the above condition, ^^^^ <^ 2 holds for 
6'sat smaller than 7r/2, for which (f ) = 0.5 < 2. 



2. Large saturation angles (small /3 and hg/l) 

For systems with small /3 the above approximation 
should fail, as is apparent from eq [32l For such cases, 
we can use a different approximation assuming that the 
gap kg is small enough, such that 

» ^ (34) 

The saturation angle can then be found from a different 
approximated form of a{^) (eg I28p. near its minimum: 



3d 



2TT£Q£dP 

3d 



1 



2/3- 



2-e 



-r 



i 

-1 



(35) 



Minimizing Q!(^), we get a quadratic equation in the vari- 
able i^: 

f -4(/3 + l)C^-f 2(/3-f 2) = 
whose solution is 

= 2(/3 + 1) ± y^4(/3 + l)2-2(/3 + 2) 



(36) 



(37) 

For small /3 it is possible to further simplify the ex- 
pression for ^sat to obtain 



3 

■sat 



2 1 




(38) 



Checking for self-consistency, the condition holds for 
small enough gaps. 
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V. SUMMARY AND OUTLOOK 

In this work we propose a novel approach towards elec- 
trowetting that, among other results, can account for 
contact angle saturation (CAS) applicable to some elec- 
trowetting setups. The model is based on a generalized 
version of the free energy accounting for various electric 
contributions. The interplay between the capillary (i^cap) 
and electric {Fd) terms depends on the applied voltage 
U, because Fd ~ —U^. Therefore, when an external 
voltage is applied it will drive the system away from its 
capillary free energy minimum and towards its electric 
free energy minimum. 

Our approach is distinctly different from other views 
of electrowetting that make use of the Young-Lippmann 
formula. In our model the electric term can exhibit a va- 
riety of dependencies on the contact angle as determined 
by the exact system geometry. Particularly, if the electric 
term Fd, has a global minimum at a certain contact angle 
^sat, then for high enough voltages this angle also min- 
imizes (asymptotically) the total free energy, Jtot- Ad- 
ditional increase of the applied voltage does not change 
the location of the global minimum and the contact angle 
saturates at 6'sat • We identify exactly this angle with the 
saturation angle found in experiments. This very general 
assumption {Fd with a minimum) is all that is needed 
to show that in the low- voltage limit a Young-Lippmann 
compatible behavior is expected, while in the high- 

voltage limit a ^ U^'^ saturation should be present. Nu- 
merical calculations suggest that combination of these 
two limiting behaviors approximates rather well the full 
expression for 0{U) in the whole voltage range, in the 
whole voltage range. 

When applying our approach to EWOD setups, we 
take two contributions to Fd into account: (i) the 
double-layer at the drop/substrate interface; and, (ii) 
another double-layer at the drop/counter-electrode inter- 
face. The latter was previously unaccounted for because 
it was considered to be negligible due to geometry, or that 
its relaxation time was considered to be very fast. How- 
ever, we estimate the relaxation time to be long (on the 
order of seconds) and, therefore, the effect of the counter- 
electrode double-layer cannot be neglected for AC sys- 
tems. Similarly, it cannot be neglected for low voltages 
(which will exclude electrochemical processes from tak- 
ing place at all) ^2l\ , or in DC applications that include a 
dielectrically coated counter-electrode. Using AC circuit 
analysis we show that Fd indeed has a global minimum 
that produces the CAS effect. 

The value of the saturation angle as well as the en- 
tire electrowetting curve 9 (U) can be found numeri- 
cally for any choice of system parameters, and our spe- 
cific choice is inspired by the experiments reviewed in 
Ref IT^ . There is a qualitative agreement with experi- 
mental results, which includes an initial compliance with 
the Young-Lippmann formula (scaled correctly) , followed 
by a cross-over to CAS. The values obtained for the sat- 
uration angle, cross-over voltage and saturation voltage 



are also compatible with experimental values. 

In addition, we investigated the frequency dependence 
of electrowetting. It is shown that the value of the sat- 
uration angle is independent of the AC frequency for a 
large range of frequencies (1 kHz to 100 kHz) for a specific 
choice of parameters. A numerical analysis of the fre- 
quency dependence of electrowetting was conducted for 
a set of system parameters inferred from RefiSj, and shows 
semi-quantitative agreement with the experiment. These 
results show that an approximation of the free energy can 
be justified, such that the entire frequency dependence is 
captured in a scaling factor of the applied voltage. It pre- 
dicts that the electrowetting effect should diminish with 
rising frequency, as indeed found in experiments 0. 

In its DC limit our model can converge to the Young- 
Lippmann formula, depending on the values of the Young 
and saturation angles. We use this result to show a novel 
way to extend the Young-Lippmann formula from DC to 
AC systems. We conclude that the validity of the Young- 
Lippmann formula is related not to the range of applied 
voltages, as it is commonly viewed, but rather to the ac- 
cessed range of contact angles. In commonly used EWOD 
setups, which are intentionally devised to have as high a 
Young angle and as low a saturation angle as possible 
(so the effect can be more easily measured), our model 
is compatible with a compliance to the Young-Lippmann 
formula at low voltages (and hence high contact angles). 
We note that the DC limit of our model can be most use- 
ful in DC applications that employ low voltages and/or 
include a dielectrically coated counter-electrode. 

Our model does not rely on any leakage mecha- 
nisms to predict CAS. Nevertheless, we would like to 
stress that leakage mechanisms treated in previous works 
[7,22,26,31-33] can be added. Interestingly, it is conceiv- 
able that a cross-over between inherent CAS (as in the 
present model) and CAS originating from leakage mech- 
anisms is responsible for the time-dependent saturation 
angle reported in Ref [U. 

The fact that the saturation angle depends on elec- 
tric parameters whereas the Young angle depends on the 
capillary parameters leads to the surprising possibility 
of reversed electrowetting. Therefore, it may be possi- 
ble to construct a system in which the Young angle is 
lower than the saturation angle. In such a system the ef- 
fect of applying an external voltage would be an increase 
in the contact angle — in total contradiction with the 
Young-Lippmann formula that allows only a decrease in 
the contact angle. We give an example of a choice of 
parameters that should yield reversed electrowetting. 

Recently, the separate control of the Young and satu- 
ration angles was demonstrated in experiments 14611. This 
ability was utilized to construct a set of dye cells [431 that 
are 'complementary' in their opposite response to applied 
voltage (black-to- white or vice- versa). We believe that 
further research in this direction will provide ample op- 
portunity to test for the existence of reversed electrowet- 
ting. Finding such evidence would have a potential im- 
pact that can go much beyond our specific model. 
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We hope that some of the predictions presented in this 
paper will be tested in future experiments in a quantita- 
tive fashion, gaining more insight on electrowetting and 
the CAS phenomenon. For example, it will be interest- 
ing to study how retracting the counter-electrode and, 
hence, reducing its contact area A2 affects the saturation 
angle, as well as coating it with a dielectric material. Our 
results suggest that more research into processes taking 
place at the counter-electrode is needed, especially with 
regards to CAS in DC EWOD setups. 



the ratio between the two capacitances can finally be ex- 
pressed as: 
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Appendix 

It is convenient to express the geometrical parameters 
and the capacitances in terms of a monotonic function of 
the contact angle 

1 a\ 1/3 

1 — cos P ^ ' 



where I — {SV/ny^^, derived from the drop volume V, 
is a characteristic length. 

The geometrical parameters defined in eq [15] can then 
be written as 

h = im 

A2 = 2nhhg {Ih-^i - 1) (A.2) 

Combining the above expressions with the definitions 
of the two capacitances, we obtain 

C2 = 2nbl^ - l-'hg) (A.3) 
With the use of a dimensionless parameter 
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